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Abstract 

We present a lattice QCD determination of the b quark mass and of the B and B s 
decay constants, performed with Nf = 2 twisted mass Wilson fermions, by simulating 
at four values of the lattice spacing. In order to study the b quark on the lattice, 
two methods are adopted in the present work, respectively based on suitable ratios 
with exactly known static limit and on the interpolation between relativistic data, 
evaluated in the charm mass region, and the static point, obtained by simulating the 
HQET on the lattice. The two methods provide results in good agreement. For the 
b quark mass in the MS scheme and for the decay constants we obtain m^imi)) = 
4.29(14) GeV, f B = 195(12) MeV, f Bs = 232(10) MeV and Jbs/Ib = 1.19(5). As 
a byproduct of the analysis we also obtain the results for the fn and fo s decay 
constants: f D = 212(8) MeV, f Ds = 248(6) MeV and f Ds /f D = 1.17(5). 
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1 Introduction 

The study of physical processes involving the b quark are of utmost importance for accurate 
tests of the Standard Model and for searching New Physics effects. On the experimental 
side, B-factories have played a fundamental role in the achievement of the present accuracy 
and further improvements are expected and looked forward from LHCb and the planned 
SuperB factories. It is therefore crucial to have theoretical uncertainties well under control, 
in particular those of the hadronic parameters computed on the lattice. 

Two particularly important cases of study are the purely leptonic decays B — > tv t and 
B s — > fi + {i~. The first process is particularly sensitive to potential New Physics contri- 
butions mediated, at tree level, by charged Higgs. The relevant entries in the Standard 
Model prediction for the decay rate are the CKM matrix element V u b, which can be ex- 
tracted from the study of semileptonic B — > n £ i/£ decays without significant New Physics 
contributions (for £ = e,fi), and the pseudoscalar decay constant fs- The measured values 
of the B — > tv t decay rate deviate, at present, by about 3 sigma from the Standard Model 
prediction [Tj [2], within relatively large experimental and theoretical uncertainties. In this 
respect, improving the lattice determination of fs would be an important ingredient for 
increasing the chances of detecting the contribution of New Physics effects to this decay. 
Another golden process for the detection of potentially large New Physics contributions is 
the rare leptonic decay B s — > ^ + fi~, which is being studied with unprecedented accuracy 
at LHCb. In this case, the relevant hadronic parameter to be determined on the lattice, 
which enters the theoretical prediction of the decay rate, is the pseudoscalar decay constant 
/b s . The determination of both f B and /# s , together with a prediction for the b quark 
mass rub, are the scope of the present study. 

With the available computer power it is not possible to simulate quark masses in the 
range of the physical b mass keeping, at the same time, finite volume and discretization 
effects under control. In order to circumvent these problems, many different methods have 
been proposed so far (see [3] for a recent review). 

In [4] we performed an exploratory calculation of the b quark mass and the decay 
constants Jb and fs s by introducing suitable ratios having an exactly known static limit. 
In [5] a more standard method [6] was applied, using lattice QCD data with the heavy 
quark mass ranging from the charm region up to more than 3 GeV, together with the 
information coming from a calculation in the static limit point. In the following, we will 
refer to the two approaches as to the "ratio method" and the "interpolation method" 
respectively. 

Here we update and finalize both the analyses, by implementing several improvements. 
We replace the preliminary values of the quark mass renormalization constants with the 
published results of j7]. We increase for some ensembles the statistics and we use more 
data, in particular, data at the finest lattice spacing (/3 = 4.2) are now included also 
in the analysis with the ratio method. The main improvement in the analysis based on 
the interpolation method consists in studying the dependence of the decay constants on 
the quark masses, instead of the meson masses, and performing the extrapolation to the 
continuum limit at fixed (reference) values of the heavy quark mass. This allows us to 
better disentangle discretization effects from the (physical) heavy quark mass dependence. 



The use of the quark masses in the determination of /# and fs s requires as input the value 
of the b quark mass, m^, which we obtain from the ratio method. 

Our results for the b quark mass in the MS scheme and for the decay constants (the 
latter obtained by averaging the results of the two methods) read 

m b (m b ) = 4.29(14) GeV , 
f B = 195(12) MeV , f Ba = 232(10) MeV , ^ = 1.19(5) , (1) 

JB 

which are in good agreement with our previous results jU [5], but which have smaller 
uncertainties. In particular, with respect to the result YnbifrTb) = 4.63(27) GeV of [3], we 
obtain for the b quark mass a central value which is smaller by approximately one standard 
deviation and a reduction of the uncertainty by almost a factor two, mainly because of the 
improvement in the determination of the quark mass renormalization constant. 

As a byproduct of the analysis we also obtain the results for the fjj and fz> s decay 
constants 

f D = 212(8) MeV , f Ds = 248(6) MeV , ^ = 1.17(5) , (2) 
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which update and improve our previous determination [8j. 

2 Simulation details 

The calculation is based on the Nf = 2 gauge field configurations generated by the Eu- 
ropean Twisted Mass (ETM) Collaboration with the tree-level improved Symanzik gauge 
action [9] and the twisted mass quark action [10] at maximal twist, discussed in detail 
in [11]- [16]. We simulated Nf = 2 mass-degenerate dynamical quarks, whose mass is 
eventually extrapolated to the physical isospin averaged mass of the up and down quarks, 
m u /d- The strange and charm quarks are quenched in the present calculation. In our lattice 
setup all physical quantities are 0(a) improved [IB] , in particular cutoff effects related to 
the heavy quark mass fih are of order a 2 fi\. 

For further details of our simulations we refer to [17] . where the same ensembles of 
gauge configurations were used. We recall here that data at four values of the lattice 
coupling, (3 = {3.80,3.90,4.05,4.20}, are included in the analysis. The corresponding 
values of the lattice spacing, a = {0.098(3), 0.085(2), 0.067(2), 0.054(1)} fm, have been 
determined in [TTJ together with m u / d using the physical values of the pion mass and decay 
constant as input. From [T7] we also take the values of the average up/down and the 
strange quark masses, namely m u /d(2 GeV) = 3.6(2) MeV and m s (2 GeV) = 95(6) MeV. 
For the quark mass renormalization constants Z^ = Zp l we use the results obtained in [7], 
i.e. Z P (MS, 2 GeV) = {0.411(12), 0.437(7), 0.477(6), 0.501(20)} at the four beta values (see 
also [T7] for the estimate of Zp at /3 = 4.20). 

At variance with [17J, where only the light, strange and charm quark masses were 
studied, a wider range of values for the valence quark masses is considered here, in order 
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Table 1: Values of simulated bare quark masses in lattice units, for the four (3 values, in 
the light (a\ii), strange (afj, s ) and heavy (apth) sectors. In the last column the minimum 
values of time t min chosen for the 2-point function fits are collected. 

to get closer to the physical b quark mass. The values of the simulated valence quark 
masses are collected in Table [TJ The values of the valence light quark mass, fii, are always 
taken identical to those of their sea counterparts. The heavy quark mass fih ranges from 
approximately m c up to 2.3-2.4m c , being m c the physical charm quark mass. Correlators 
at higher fih values have been simulated and were included in [5] . They are characterized by 
large fluctuations in the effective mass plateaux, and thus by large statistical uncertainties. 
As a consequence, these data turn out to be irrelevant in the fits, and we have excluded 
them from the present analysis. 

We now proceed to describe the two approaches adopted in the present work to study 
the .B-physics observables, namely the ratio method of [4] and the interpolation method. 

3 Ratio method 



3.1 The b quark mass 

The b quark mass is obtained by implementing the ratio method of [4J, briefly recalled 
hereafter. The method is suggested by the HQET asymptotic behavior of the heavy-light 
meson mass Mm, 
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where u.J i oe is the pole quark mass and the limit ([3]) is approached without corrections 
of 0(l/log(/^ oe /AQ CD )). The first step is to consider an appropriate sequence of heavy 
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and ranging from the charm mass to values somewhat below the bottom mass. Here and 
in the following we denote by a "bar" the quark masses renormalized in the MS scheme 
and, if not otherwise specified, a renormalization scale of 2 GeV is implied. 

Then one computes the following ratios that have an exactly known static limit, 
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The function p(/x/i, //*) is the factor that converts the renormalized MS quark mass (at the 
scale /i*) into the pole mass, 
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p h = p(Ph,p )Ph{p ) , 

known up to N 3 LO in perturbation theory [HI [19]. The NLO expression reads 
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used for Nf — 2 in the present analysis. We notice that the dependence on the scale \i* 
cancels in the ratios of p factors evaluated at different heavy quark masses and thus in the 
y ratio defined in eq. (jSJ). 

From eq. ([3]) and QCD asymptotic freedom it follows that the ratios (jSJ) have the 
following static limit: 

lim y(p h , A; fy, a = 0) = 1 . (8) 

The value of the ratio A of eq. (jlj), between two subsequent values of the heavy quark 
mass, is chosen in such a way that after a finite number of steps the heavy-light meson 
mass assumes the experimental value Mb = 5.279 GeV (we find A = 1.1762) . In order to 
implement this condition, the lattice data at the four lattice spacings are interpolated at 
the following values of the heavy quark mass, 

1.140 GeV , pf = \pjjp = 1.341 GeV , pf = X 2 p^ = 1.577 GeV, (9) 

1.855 GeV , u[. 5) = A 4 /!^ = 2.182 GeV , u ( , 6) = A 5 /!^ = 2.566 GeV. 



rf' 



-(4) 
Ph 



A p h 



Ratios of the kind defined in eq. ([5]) are introduced because, besides having an exactly 
known static limit, they are also expected [4] to have a smooth chiral and continuum limit, 
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1: Light quark mass dependence of the meson mass Mht(fA ) (left) and of the ratio 
(right) at the four values of the lattice spacing. 



as shown in the right plot of fig. [TJ In fig. [JJ (left) we show the chiral and continuum 
extrapolation of the heavy-light meson mass evaluated at the first of our heavy quark 
masses pr h , namely Mhi(p> h ), at the four available f3 values. We have considered the 
following (phenomenological) linear fit which, as shown in fig. [TJ turns out to describe well 
the lattice data, 



M hl (jx 
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After performing the continuum and chiral extrapolation of the ratios (JS)), we study 
their dependence on the inverse heavy quark mass. Inspired by HQET, we perform a 
polynomial fit in 1/jS/i, of the form 



y(ji h ) = 1 + — + — , 
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which imposes the constraint y — 1 at the static point. The fit is illustrated in fig. [2 
A detailed discussion of the \ih dependence of the ratio y, in comparison to the HQET 
expectation, is provided in the Appendix. 

The value of the b quark mass is finally determined by considering the following equation 
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which should be looked at as a relation between the mass of the heavy-light meson, 
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and the corresponding heavy quark mass fi h , being Mh u /d(p 
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initial triggering value. The b quark mass is then determined by finding the value of K at 
which Mhu/diflh ) takes the experimental value of the 5-meson mass, Mb- Calling K^ 
the solution of the resulting eq. ( 112"]) (we find Kj, — 9), one gets for fit, = mb(2 GeV) the 
simple relation 

jj b = A^4 1} = 4.91(15) GeV . (13) 



1.01 




0.00 0.10 0.20 



0.30 0.40 0.50 
1/W, (GeV- 1 ) 



0.60 0.70 0.80 



Figure 2: Heavy quark mass dependence of the ratio y(fth) extrapolated to the physical 
value of the light quark mass and to the continuum limit. The vertical line represents the 
value of the physical b quark mass. 



We observe that it is always possible to guarantee that the solution Kb is an integer number 
through a slight variation of the parameter A and/or of the triggering mass p, h . 

An equivalent method consists in determining the b quark mass by studying Mh s instead 
of Mh u /d and using in input the experimental 5 s -meson mass value, Mb s — 5.366 GeV. A 
very similar result is obtained from this analysis: \x\> = 4.92(13) GeV. The small difference 
(0.01 GeV) with respect to eq. ( TT31 indicates a good control of the chiral extrapolation 
which, in particular, in the heavy-strange meson case involves only the sea quark mass. 
The main effect of the uncertainty due the chiral extrapolation is accounted for by the 
error quoted in eq. ffT3|) . which comes from the chiral, continuum and l/fih fits. 

In order to estimate the residual uncertainty due to discretization effects, we have 
tried to include in the continuum extrapolation, besides the leading 0(a 2 ) correction, an 
additional a 4 term. We find, however, that this subleading contribution cannot be fitted 
with our data. Therefore, we have repeated the analysis by excluding the data at the 
coarsest lattice spacing (/3 = 3.80). The difference in the determination of the b quark 
mass turns out to be of 0.05 GeV. 

In order to estimate the systematic error associated to the interpolation of y{jih) a s a 
function of 1/fih, we have repeated the whole analysis by choosing a third order polynomial 
in 1/fih (rather than a second order one, as in the ansatz (1111) ). This change resulted in 
an increase of the b quark mass of about 0.5%, corresponding to a shift of ~ 0.02 GeV of 
the central value result of eq. ffTB"]) . 

An additional uncertainty is introduced by the truncation of the perturbative series in 
the determination of the pole mass in eq. ([6]), which is affected by renormalon ambiguities. 
When comparing the results obtained with the NLO definition of the pole mass to the 
results found with the LO one, the difference in the b quark mass is found to be small, 
of about 0.01 GeV. The sensitivity to the pole mass definition, which appears in the 
intermediate steps, thus largely cancels out in the final determination. 



We finally quote the b quark mass at the conventional renormalization scale of m b itself 

m b (m b ) = 4.29(13) (4) GeV =4.29(14) GeV, (14) 

where the first error is of statistical and fitting origin and the second one is the sum in 
quadrature of the residual systematic uncertainties discussed above. The present result 
for the b quark mass has a central value which is smaller by approximately one standard 
deviation than the value found in [4], and an uncertainty which is reduced by almost a 
factor two, reflecting the various improvements implemented in the present analysis. 

3.2 Decay constants 

A similar strategy is employed by applying the ratio method to determine the B and B s 
meson decay constants. The HQET asymptotic prediction for the decay constant is 
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Therefore, the ratios with static limit equal to one of interest in this case are, for fs and 

,- x - ) = x i/2 fhi{jih,iii,a) Ci at (nl,ji h /X) [p(fih, /j*)] 1/2 
z^ h ,A,Pt,a) _ fafa/x^a)' Cr^lH) [p(/i,/A,^)]V2 



z s{n-hi A; Hi, fa, a ) 



A i/ 2 fhs(Vh,Pe,Ps,a) Cj at (p* b ,p h /X) [p{p h , /i*)] 1 / 2 
fh s (fih/\P>t,fis,a) C^ifilfih) [pifih/X,/!*)} 1 / 2 ' 



The ratio of p factors (raised to the appropriate power) is present to convert MS heavy 
quark masses to pole masses as in eq. (jSJ). The factor C^ at (pl,flh), defined as 
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provides the matching between the decay constant in QCD for a heavy quark mass fih and 
in HQET, and the running of the static axial current to the renormalization scale pi, and 
it is known up to N 2 LO in PT [20]. The NLO expression used in the present analysis reads 
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In order to have better control on the chiral extrapolation, we consider as primary 
quantities in the present analysis the decay constant fs s , whose dependence on the light 
quark mass only occurs through sea effects, and the ratio /b s //b which provides a direct 
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Figure 3: Light quark mass dependence of the decay constant fhaifA ) (left) an d of the 
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ratio z s (p h ) (right) at the four values of the lattice spacing. 
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Figure 4: Light quark mass dependence of the ratio of decay constants fhs(lAi )l fhi{lA 
(left) and of the double ratio z a (pr h )/z(pr h ) (right) at the four values of the lattice spacing. 



determination of the SU(3) breaking effect in the decay constant. Within the ratio method, 
these quantities are obtained from the ratio z s and the double ratio z s /z. 

Both z s and z s /z have a smooth chiral and continuum limit, as illustrated in figs. |3] 
and|H In particular, the results for the double ratio z s /z turn out to be well described by 
both a linear and a constant behavior in both fy and a 2 (see fig. H] right). For simplicity 
reasons the constant fit ansatz was chosen. In the left panels of figs. [3] and H] we show the 
chiral and continuum extrapolation of fh s and fhs/fhi at the initial (triggering) mass pr h . 
F° r fhs/ fhi, heavy meson chiral perturbation theory (HMChPT) predicts at the NLO a 
linear+logarithmic dependence on the light quark mass, since a chiral log controls the chiral 
behavior of fs (see eq. (j20|) below). With our results, the logarithmic dependence cannot 
be appreciated, and we thus perform also a simpler linear fit in the light quark mass which 
turns out to describe well the lattice data. As discussed in section 5, we eventually average 
the results obtained from the HMChPT and the linear fits and include the difference in 
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Figure 5: Heavy quark mass dependence of the ratio z s (fih) (left) and of the double ratio 
z s {fih)/z{fih) (right) extrapolated to the physical value of the light and strange quark 
masses and to the continuum limit. The vertical line represents the value of the physical 
b quark mass. 



the systematic uncertainty. For fh s , which depends on the light quark mass for sea effects 
only, we have implemented both a linear and a quadratic fit. They turn out to provide 
essentially identical results. 

Finally, we study the dependence of the ratio z s and the double ratio z s /z on the heavy 
quark mass, which are shown in fig. |5j For z s we perform a quadratic interpolation to the 
h quark mass as for the ratio y{fih) in eq. ( flTT) and, also in this case, a detailed discussion 
is provided in the Appendix. For z s /z, the dependence on the heavy quark mass is barely 
visible, so that in this case we perform either a linear interpolation in 1/fih or we fix this 
ratio equal to its asymptotic heavy-quark mass limit, z s /z = 1. 



4 Interpolation method 

As already mentioned, the interpolation method consists in interpolating to the b quark 
mass the relativistic results obtained for values of the heavy quark masses in the range 
around and above the physical charm (up to twice to three times its value) and the result 
evaluated in the static limit by simulating the HQET on the lattice. In this section, 
we describe these results by addressing, in turn, the calculation with relativistic lattice 
QCD in the charm mass region, the calculation within the HQET on the lattice, and the 
interpolation among the two sets of results. 



4.1 Decay constants in relativistic QCD 

The lattice relativistic data for the heavy-light and heavy-strange meson masses and decay 
constants are the same used for the ratio method. We considered in the analysis four values 
of the lattice spacing and the values of valence quark masses collected in Table [TJ With 
respect to the preliminary results with this method presented in [5], we added an ensemble 



with a lighter quark mass at /3 — 4.2 and, for other ensembles, we increased the statistics. 
Another update w.r.t to the analysis in [5] concerns the renormalization constants, which 
had preliminary values at the time of [5], and have been later updated and published 
in [7]. The main improvement, however, concerns the disentanglement of the heavy mass 
dependence from discretization effects. In the present analysis the extrapolation to the 
continuum limit is performed at fixed (renormalized) heavy quark mass. The whole analysis 
consists in the following steps. 

First, we slightly interpolate the lattice data to reach a set of reference heavy quark 
masses equal at the four /3 values. This allows us to study discretization effects at 
fixed heavy quark mass. We have chosen the reference masses within the range of the 
simulated values. In the MS scheme, at ji — 2 GeV, the set of reference masses is 
pL h = {1.25, 1.50, 1.75, 2.00, 2.25} GeV. 

We have then performed a combined continuum and chiral extrapolation, at fixed ref- 
erence heavy quark mass. As for the ratio method, we consider as primary quantities fh s 
and the ratio fhs/ f he obtained in the present analysis from & hs and &hs/&he, where 



®h s = fhsVMhZ and Zh±= J -±± x ^±. (19) 



ffts^A, Mh 

$ M f he V M h 

An alternative analysis based on the definition of &hi(s) i n terms of the pole mass, rather 

than the meson mass, i.e. $«( s ) = fhe(s)\I Vft° e > nas Deen also performed, leading to fully 
equivalent results. 

The light quark mass dependence predicted for $^ and &hs by HMChPT [2] 
at the NLO reads 

* ( - - \ A U 3(1 + 3^ 2 ) 2B fi £ f 2B fl t \ _ 2 

^ hs (a,fj, E ,fx s ,fJL h ) = D h (1 + E h Jx e + F h fx s + G h a 2 ) , (20) 

where we have also included in the above expressions a linear dependence on a 2 to account 
for leading discretization effects. The subscript h in the fit parameters of eq. ( l20j) denotes 
the dependence on the heavy quark mass. 

As previously discussed for the ratio method, the contribution of chiral logs in the ratio 
^hs/^Mi predicted by HMChPT, cannot be appreciated with our data (see fig. H] left). 
Thus, in order to perform the chiral extrapolation to the physical light quark mass, we also 
perform a simple linear fit in p,£ and eventually take the average of the two results. In the 
fit based on HMChPT, we take for the parameter g the value g = 0.61(7) obtained from 
the experimental measurement of the gD*Dn coupling [21]. We choose this value, instead of 
the HQET prediction g = 0.44(8) [251126] , as we fit data that are close to the charm mass 
region and in order to conservatively include in the average the maximum spread resulting 
from the different ways of performing the chiral extrapolation of our data. For $^ s , as for 
fhs within the ratio method (see fig. [3] left), we have tried both a linear and a quadratic fit 
in ji£, obtaining very similar results. 

For illustration, the size of discretization effects in the calculation of $^ s is shown in 
fig. El for a simulated value of the light quark mass of about 50 MeV and with fi s rj Ji p s hys . 
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Figure 6: Dependence of §hs on the squared lattice spacing, at the five reference values of 
the heavy quark mass (for a fixed value of the up/down and strange quark masses). 

It is interesting to note that lattice artifacts turn out to be small. We find that this is a 
consequence of a partial cancellation between discretization terms in the decay constant 
and in the rooted meson mass, which are of similar size and opposite sign. For the same 
reason, the ratio $hs/&hi turns out to be practically independent of the lattice spacing. 

4.2 Decay constants in the static limit of HQET 

We now summarize the procedure adopted to evaluate the static-light meson masses and 
decay constants, which follows and improves the analysis of [5]. We use the same setup as 
for the recent ETMC study of the static-light meson spectrum [27J [2S]. Technical details 
regarding e.g. number of ETMC Nf = 2 gauge configurations considered, meson creation 
operators, corresponding correlation matrices and their analysis, smearing techniques to 
enhance the signal quality [29] as well as efficient propagator computation can be found in 
these references. 

The lattice action used to describe the static quark is the HYP2 static action [2"9"| I3"0"| I3"T] 

S h = a 4 ^2i> h (x)D ip h (x) = a 3 ^2$ h (x)[ip h (x) - V^ YP (x - a6,0)ip h {x - a6)] , (21) 



where Vhyp is the so-called HYP-link, which is a gauge covariant function of the gauge 
links located within a hypercube. 

Static-light correlators have been calculated for a subset of the configuration ensembles 
used for the relativistic calculation, namely for two values of the lattice spacing, /3 = 3.90 
and = 4.05. 

Due to the mixing between the pseudoscalar and the scalar currents occurring in the 
static-light framework, $§/*•, is obtained as a linear combination of two matrix elements, 



nf = Z^(0\A l5 x q \Bff) +ir q Z^(0\A Xq \Bff) 



(22) 
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Figure 7: Tie combination $gf /Z(4.05) (2eft), with Z = (Z P tat + Zf tat ) /2, and the ratio 

(TjStat 

^|t (right) as a function of the (squared) light pseudoscalar meson mass. In the left plot 
the data at (5 = 3.90 have been multiplied by the appropriate factor to match the same 
scale for the data at ft = 4.05. 



where Xq = e~ %rq ^^^ 5 il) q is the light quark held in the twisted basis, r q = ±1 is the corre- 
sponding Wilson parameter, and Zp at and Z| tat the static-light renormalization constants. 
In order to improve the signal quality of the lattice data we have combined the results with 
r q = +1 and r q = — 1, which are identical on average due to twisted mass symmetries. 

In order to check the stability of the fitting method to extract the above matrix elements, 
we have performed, as in [27J, the computations with different parameters: number of states 
to be fitted in correlators, fitting time ranges, operator content of the correlation matrices. 
We have obtained results that are consistent within statistical errors. 

For our setup (tree-level improved Symanzik gauge action and HYP2 static action) 
the non-perturbative values of the static-light renormalization constants are not available 
at present. We thus rely on the perturbative estimate at one loop [32J evaluated using 
a boosted coupling. It provides Z P tat =J0.85(8), 0.86(7)} and Z| tat = {0.93(4), 0.94(3)}, 
at (3 = {3.90,4.05} respectively, in the MS scheme at the renormalization scale 1/a. The 
uncertainty on Zp*g has been conservatively estimated as half of the deviation from unity. 

The chiral extrapolation of $§ at and cDgaty^&tat j^ i oeen performed, as for the rela- 
tivists data, using the HMChPT functional forms given in eq. (120 p . and it is shown in 
fig. [7] in terms of the light meson mass squared Mf t . In this case, however, discretization 
terms are set to zero, as the fit includes data at only two values of the lattice spacing. 

The matrix element shown in the left plot of fig. [7] is chosen to be, for better clarity, 
the ratio $§ at /Z, where Z = (Z| tat + Zp tat )/2. This quantity only depends on the ratio of 
renormalization constants, Z| tat /Zp tat , but not on their average. One can also show that 
the relative systematic errors due to inaccurate knowledge of Z| tat /Zp tat on these matrix 



elements are A(Z| tat /Z|> tat ) x 0(a), where A(Z| tat /Z 



stat \ 



is the error on Z| tat 



/Z 



p at and 



O{o) denotes matrix elements which are proportional to the lattice spacing. Consequently, 
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the matrix elements shown in fig. [7J are only slightly affected by the systematic uncertainty 
on the renormalization constants due to the use of one-loop lattice perturbation theory. Of 
course the evaluation of $|£f at each lattice spacing requires the info on one renormalization 
constant, which is conveniently chosen as (Z| tat + Zp at )/2. 

Fig. [7] shows that the chiral dependence on the light pseudoscalar mass is found to be 
very smooth for ^f^f, and a simple constant fit is used to perform the chiral extrapolation. 
In order to estimate the uncertainty due to the chiral extrapolation, a linear fit is tried for 
comparison and the difference is included in the systematic error. For the ratio cDgay^gat^ 
besides the fit based on HMChPT we have also considered both a linear (in Mf h ) and a 
quadratic fit. The variation of the result is found to be of the order of 5% and it is included 
in the systematic uncertainty. 

It is clear from fig. [7J (right) that the data for the ratio (jj^y^stat are a ff ec ^ed by 
negligible cut-off effects. Some discretization effects are visible, instead, in the left plot of 
fig. [7J which shows the results for $§if. In order to estimate their magnitude we compare 
the results for $f£f obtained by either fitting together data at /3 = 3.90 and /3 = 4.05 or 
fitting data at one lattice spacing only. The difference, which turns out to be at the level 
of 3%, is included in the systematic uncertainty. 

For $§ at the simulated (bare) strange quark mass is fixed to the values a/j, s = 0.0220 at 
[5 = 3.9 and a/i s = 0.0170 at (3 = 4.05 respectively. These values correspond to the physical 
strange mass, as obtained from the analysis at fixed lattice spacing. The continuum limit 
performed in the study of [UJ has later provided a value for the strange quark mass which 
is smaller by approximately 22% at /3 = 3.9 and 13% at /3 = 4.05. The effect of this 
mismatch is discussed in the following section and included in the systematic uncertainty. 

As a last step in the analysis of the static data, a perturbative evolution at NLO [20] 



has been applied to evolve the results in the MS scheme from the initial scale \i = 1/a to 
a common reference scale /i£ = 4.5 GeV, obtaining 

&g?(l4) = 0.67(4) GeV 3/2 and $g a y $ stat = h2 s(7) . (23) 

For the B meson, eq. ( 1231) corresponds to the result 

$%*(£) = 0.52(3) GeV 3/2 . (24) 

As a further consistency check between the analyses based on the ratio and the interpo- 
lation methods, we have used the ratio method to predict the value of $b s /$b in the static 
limit, finding $g^y$g at = 1.20(5). The latter is compatible with the result in eq. (|23|) 
obtained from the direct lattice simulation in the HQET. 

4.3 Interpolation of relativistic and static data 

In order to perform a combined fit of relativistic and static data, we convert the relativistic 
values of $>hs from QCD to HQET, by using the NLO matching and evolution factor C^ at 
(see eqs.( lT7|) and (TT8|) ). The renormalization scale is chosen to be [i\ = 4.5 GeV as for the 
static data. Note that in the ratio Qhs/^M the C^ at factor cancels out. 

The interpolation is then performed, as shown in fig. [HI through a fit in 1/fih, which 
is quadratic for §hs (similarly to eq. (TTTT) ) and only linear for $/i S /$w, where a much 
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Figure 8: Dependence of$>hs (left) and <&hs/®M (right), m the chiral and continuum limit, 
on the inverse of the heavy quark mass. 

smoother dependence on the heavy quark mass is found, as expected. Finally, the physical 
results for the decay constants are obtained by inserting the "physical" value of the b quark 
mass determined from the ratio method, given in eqs. (113j) - (114p . The uncertainty on the b 
quark mass has been propagated by assuming a gaussian distribution of the errors. 

5 Results for the decay constants 

In this section we present and discuss the final results obtained for the decay constants 
fss, Ibs/Ib and f B from the ratio and the interpolation methods. 

As discussed in the previous sections, in order to estimate the uncertainty due to the 
chiral extrapolation we compare the results of two different chiral fits. This comparison is 
relevant in particular for the ratio <&hs/<&he, from which &bs/&b is extracted. In this case, 
the two fits are based either on the linear + logarithmic dependence on the light quark 
mass predicted by HMChPT (eq. (1201) ) or on a polynomial (quadratic) behavior. For $ hs 
we have tried both the linear fit of eq. ( 120]) and a quadratic fit. The analogous chiral fit 
ansatz employed in the analysis of fhaifihj/fheifih) within the ratio method framework 
have been discussed in section 3.2. The results are collected in Table |2] for both the ratio 
and the interpolation method. The first error quoted in the table is the one coming from 
the fit, and includes both the statistical error and the systematic uncertainty due to the 
chiral and continuum extrapolation and to the interpolation to the b quark mass. The 
second error accounts for the additional systematic uncertainties and it has been evaluated 
as follows: 

- Continuum limit : when performing the continuum limit, both in the ratio and the 
interpolation method, we consider a linear fit in a 2 . Since an additional a 4 term 
cannot be fitted with our data, we estimate the uncertainty due to discretization 
effects by excluding data at the coarsest lattice (ft = 3.80). The central values for 
fss change by 2 and 1 MeV for the ratio and the interpolation method respectively. 
The corresponding changes in the values of the ratio JbsIJb are instead negligible. 
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f B .[ MeV] 


Ibs/ ' Jb 


Ratio Method 


Interpol. Method 


Ratio Method 


Interpol. Method 


Lin. 


Quad. 


Lin. 


Quad. 


HMChPT 


Polyn. 


HMChPT 


Polyn. 


225(7) (4) 


225(7) (4) 


237(9) (4) 


238(9) (4) 


1.22(2)(0) 


1.14(2)(0) 


1.22(5)(2) 


1.16(6)(2) 


225(7) (4) 


238(9) (4) 


1.18(2)(4) 


1.19(5)(3) 


232(10) 


1.19(5) 



Table 2: Collection of the results obtained for fg s and fss/fB from the ratio and interpo- 
lation methods. The statistical and systematic uncertainties are summed in quadrature. 
The third and fourth lines provide info on the results obtained by extrapolating to the 
physical pion mass point by using different chiral fit ansatz (see text). The final values, 
given in the last row, are an average of the results of the two methods. 

- Heavy mass dependence: Within the interpolation method we estimate the uncer- 
tainty in reaching the physical bottom mass by including, for each fl, data at two 
larger values of fih, and by choosing slightly different values for the reference masses. 
We find that with these variations the central values obtained for /# s change by ap- 
proximately 3 MeV, while the results for the ratio JbsIJb are practically unaffected. 
In the context of the ratio method analysis in order to estimate the systematic error 
associated to the determination of z s {ji^) we have varied the fit ansatz by considering 
either a second order or a third order polynomial in 1/fih- This change produces only 
a 1 MeV decrease in the final value of Jb s — 225 MeV (see Table |2]). Even smaller 
is the relative uncertainty in the 1/flh interpolation of the double ratio z s (p,h)/z(p,h), 
owing to the very flat profile of data within errors, as it is seen from fig. [5] (right). 

- Pole mass : as the pole mass is affected by renormalon ambiguities, in the analysis 
based on the ratio method we compare the results obtained by using the NLO def- 
inition of the pole mass to the results found with the LO definition. Within the 
interpolation method, instead, we have also considered the alternative definition of 
$he(s) m terms of the pole mass (rather than the meson mass), again using either 
the NLO or the LO definition of the pole mass. In both cases, the differences are 
found to be small, at the level of 1 MeV, for the decay constants, as the sensitivity to 
the pole mass definition, which appears in the intermediate steps of the calculation, 
largely cancels out in the final determinations. The results for the ratio JbsIJb are 
practically unaffected. 

- Mismatch of the strange quark mass in the static simulation: as discussed in section 
4.2, the static-strange correlators have been calculated with a value of the strange 
quark mass that was estimated from an analysis at fixed lattice spacing, and turned 
out to be larger with respect to the continuum limit estimate by approximately 22% 
at (3 = 3.9 and 13% at (3 = 4.05. In order to evaluate the systematic uncertainty due 
to this mismatch, we have analyzed the relativistic data for $^ s which are available 
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at several values of the strange quark mass. By using the continuum estimate of the 
strange quark mass, $s s decreases by approximately 2%. A similar effect can be 
thus expected for the static data. We have thus repeated the interpolation to the b 
quark mass using for the static points results smaller by 2%. We find that the B s 
decay constant decreases by 3 MeV and /s s //s by 0.015. We conservatively ignore 
the sign of the variation and consider these changes as a symmetric contribution to 
the systematic uncertainty. This uncertainty does not affect the ratio method, since 
in this case the static limit of z s and z s /z is exactly known. 

For both methods we add in quadrature the systematic uncertainties and, finally, as 
shown in the last row of Table [2J we average the results of the two methods obtaining 

f Bs = 232(10) MeV, ^ = 1.19(5) , (25) 

JB 

and for the B decay constant, which is determined for each analysis as fn = /bs/(/bs//b), 

f B = 195(12) MeV. (26) 

These values are in agreement and improve the results obtained in [I] and [5]. 

As a byproduct of the analysis we also obtain the decay constants for the D and D s 
mesons. In order to determine these quantities, we only consider three values for the heavy 
quark reference masses around the physical charm quark mass. By interpolating to the 
physical value m c (m c ) = 1.28(4) GeV, obtained in [5], we find 

f D = 212(8) MeV , f Ds = 248(6) MeV , ^ = 1.17(5) , (27) 

JD 

to be compared with the results f D = 197(9) MeV, f Ds = 244(8) MeV and f Ds /f D = 
1.24(3) of [8j. With respect to [8], the present analysis is improved essentially for the reasons 
discussed for fs and fss, namely: the statistics is increased for some ensembles, data at the 
finest lattice spacing (/3 = 4.2) are now included, the continuum extrapolation is performed 
at fixed (reference) heavy quark masses. Moreover, as discussed for JBs/fB, we perform the 
chiral extrapolation of fhs/ fhi either following HMChPT or a linear dependence on fi£. In 
[SJ the value Jds/ Id — 1.24(3) was obtained from the HMChPT fit only, while the result 
given in eq. $27]) is an average of f Ds /f D = 1.21(2) from HMChPT and f Ds /f D = 1.12(2) 
from the linear fit. By considering both results we have increased the uncertainty associated 
to the chiral extrapolation. 

6 Conclusions 

We have presented a lattice determination of the b quark mass and of the B and B s decay 
constants, obtained with Nf = 2 twisted mass Wilson fermions. Two methods have been 
employed, following and improving our previous analyses in [4] and [5]. 

The first method is based on suitable ratios with exactly known static limit and smooth 
chiral and continuum limit. With respect to [4], the present analysis includes data at four 
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values for the lattice spacing, a larger statistics, and uses the published values for the 
quark mass renormalization constants [7] and for the physical up/down and strange quark 
masses [P7] . 

The second method consists in interpolating between relativistic and static data. With 
respect to [5], we added one ensemble at j5 — 4.2, increased the statistics and, again, 
used the published values for the renormalization constants and light quark masses. A 
further improvement has been achieved by studying separately discretization effects and 
the (physical) dependence on the heavy quark mass. This has been done by performing 
the continuum extrapolation at fixed reference heavy quark mass. 

The systematic uncertainties due to the chiral and continuum extrapolation and to the 
interpolation to the physical b quark mass, as well as the sensitivity to the pole mass defi- 
nition, have been carefully studied. An important uncertainty affecting the determination 
of the ratio JbsIJb an d, in turn, of fs, is introduced by the chiral extrapolation to the 
physical value of the average up/down quark mass. We note, in this respect, that given 
an assumption for the chiral extrapolation fitting function, i.e. either including or not the 
leading chiral logarithm, the results obtained for the ratio /bs/Ib by using the ratio and 
the interpolation method are in perfect agreement within each other (see Table |2]). In order 
to reduce the uncertainty due to the chiral extrapolation, simulations at smaller values for 
the light quark masses, closer to their physical values, are needed. 

The difference between the results obtained for /# s by using the ratio and the interpo- 
lation method (approximately 5%, see Table |2]) provides an indication of the uncertainty 
due to the interpolation to the heavy b quark mass. In this respect, the main advantage 
of the ratio method is that the static limit of the ratios is exactly known (by definition), 
so that the approach does not require a dedicated lattice simulation within the HQET. 

The final results for the b quark mass in the MS scheme and for the decay constants 
read 

m b (m b ) = 4.29(14) GeV , 
f B = 195(12) MeV, f Ba = 232(10) MeV , ^ = 1.19(5) . (28) 

JB 

As a byproduct of the analysis we also obtain the results for the fo and fo s decay constants 
f D = 212(8) MeV , f Ds = 248(6) MeV , ^ = 1.17(5) , (29) 

JD 

which update and improve our previous determination [8]. 

The only systematic uncertainty which is not accounted for by our results is the one 
stemming from the missing strange and charm quark vacuum polarization effects. A com- 
parison of our Nf = 2 result for the B and B s decay constants, to existing results from 
Nf = 2 + 1 quark flavor simulations (33J, [33] suggests that the error due to the partial 
quenching of the strange quark is smaller at present than other systematic uncertainties. 
In this respect we mention that simulations with Nf = 2 + 1 + 1 dynamical flavors are 
already being performed by ETMC and preliminary results for several flavor physics ob- 
servables have been recently presented [351 EE] . 
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Appendix: Phenomenological analysis of the ratios y 

and z s 

The lattice results for the ratios y and z s derived in section 3 deviate from the their static 
limit value, in the whole range of heavy quark masses from the charm mass value up to 
infinity, by only a small amount. Specifically, with the chosen value A ~ 1.18, the deviation 
is not more than 1.5% and 4% for y and z s respectively. By looking at the best fit curves of 
y and z s as functions of 1/fth, see figs. |2] and |5] (left), one notices however a clear curvature, 
thus signalling a large 1/fil contribution in the heavy quark expansion compared to the 
linear term. In this appendix, we wish to show that this behavior is actually in good 
agreement with the predictions of the heavy quark expansion of Mm and fh s once one 
employs phenomenological or lattice based estimates for the relevant coefficients. 

We first discuss the phenomenological analysis of the ratio y(fth) defined in eq. ([5]). 
By introducing in the expression for y the heavy quark expansion for the heavy-light 
pseudoscalar meson mass, 

one finds 



A !^fl + f< A - + 3A V* + 1) + a 2 H >£.- ' 



V = 1 - A^ + i*™S2 (A p- + 1) + A^* ±_^i , (31) 



where A pole is a smoothly varying function of \Xh defined as A pole = ^° c ([ih)/ p^ e (Ph/ty = 

*P(fih)/p(fih/X)- 

In order to estimate the ratio y we considered the following phenomenological values 
for the HQET parameters 

A = 0.39(11) GeV, Ai = -0.19(10) GeV 2 , A 2 = 0.12(2) GeV 2 . (32) 
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Figure 9: The lattice QCD results for y (blue points) as a function ofl/jXh are compared to 
the phenomenological estimate. For the latter, the black solid and dashed curves represent 
the mean value and the one-standard deviation band respectively. 



Notice that while A2 is rather precisely determined from the B-meson mass splitting M^» — 
Mg, the values of A and Ai are inferred from the analysis of the inclusive semileptonic B- 
decays carried out in ref. [3?]. Using eq. fl32|) we obtain the phenomenological estimate of 
y shown in fig. which is compared to the lattice QCD results of section 3. 

The above analysis shows that the lattice results for y are consistent with the phe- 
nomenological estimate. In particular, one sees that the curvature of y as a function of 
l//2/j,, observed in the lattice QCD data, is due to the fact that the coefficients of the linear 
and quadratic terms of the heavy quark expansion in eq. (I3TI) are comparable in size and 
opposite in sign. As a result, the ratio y has a minimum in the region around the charm 
quark mass, which is the one covered by the lattice data. 

A similar analysis can be also carried out for the ratio z s , defined in eq. (TT6l) . In 
this case, the analysis is based on the heavy quark expansion for the pseudoscalar decay 
constant, 



®hs(fJ-h, Hi) 



(fhsVWs 



iQCD 



*o(k; 



x + gjM) + *2(P 



pole \ 2 



+ 



(z^ 3 



, (33) 



(~~istat(r, i,*\ " v " u/ 1 pole ' / pc 

where u* b is the renormalization scale in the HQET. Using this expansion (and omitting 
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for better clarity in the following the dependence on /i£), one finds 



1/2 $hs{^h 



\pole I \polc _ I 

v h \A ) 



i - ^~^~ " (^(^ polc + i) - ^ pole ) yi^r ■ ( 34 ) 



The expansion for z s can be then obtained by combining the above expression with the 
heavy quark expansion for the ratio y s . The latter has the same form of eq. (|3T|) but with 
HQET parameters A s , Ai s and \2 S depending on the strange quark mass. For the purpose 
of the present exercise we take 

A s = A + M Bs - M B , Ai, = Ai , \ 2s = A 2 . (35) 

We are not aware of phenomenological estimates of the HQET parameters $o 5 ^i and 
$2- Therefore we consider in this case a set of values inferred from the lattice results for the 
heavy-light meson decay constants obtained by the HPQCD collaboration and presented 
in [38]. Their result for the static parameter is $o — 0.60(4) GeV 3 ' 2 , which is consistent 
with our calculation within lattice HQET, $ — 0.67(4) GeV 3 ' 2 , see eq. (|23|) . The estimates 
for the parameters $i and <J>2 can be derived, in turn, by requiring that eq. (133]) provides 
the HPQCD determinations for fjj s and f Ba at the physical charm and bottom quark 
masses, i.e. f Ds = 249(2) MeV and f Bs = 224(4) MeV [38J. For this determination we also 
used the experimental values of the D s and B s meson masses [24] and the values of the 
charm and bottom quark masses obtained in section 3.1, namely fi c = 1.14(4) GeV and 
fib = 4.91(15) GeV. In conclusion we considered the values 

$o = 0.60 GeV 3/2 , $i = -0.48 GeV , $ 2 = 0.08 GeV 2 . (36) 

The resulting phenomenological estimate of z s corresponding to the above set of values 
for the HQET parameters is shown in fig. [TO] and compared to our determination from 
section 3. By also considering that a quantitative estimate of the uncertainties on the 
HQET parameters is beyond the scope of the present exercise we conclude, again, that the 
shape of the phenomenological curve, including its curvature, is well consistent with the 
lattice data. 
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